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(N 

ox): 

^ ' The one-dimensional problem of the nonlinear heat equation is considered. We 

assume that the heat flow in the origin of coordinates is the power function of time 
and the initial temperature is zero. Approximate solutions of the problem are given. 
Convergence of approximate solutions is discussed. 

> 1 

1 The problem statement. The nonlinear heat equation has the form [1] 

on ■ du d ( „du\ 

o: «=**(" r>0 - i>0 (1) 

where u(r,t) is temperature in the point r; t is time; x is coefficient of the heat 
i-£h ■ conductivity. 



Suppose the boundary condition has the form 

S U dr 

> 



-q t k , t>0 (2) 



=o 



^ ■ of coordinates. 



Where go is constant. Equation (J2J) corresponds to the energy flow in the origin 



We also take the initial condition in the form 

u(r,t = 0)=0, r>0 (3) 

The boundary value problem (0) — (J3J) was considered [1-4] . When the tempera- 
ture specified on the boundary the problem for equation ([I]) were considered in [3-7] . 
Numerical solution of the problem ((TJ — can be found using difference method [8] . 

In this work we are going to look for approximate solutions of the problem 

©• 

Using 



U = V n 



we get 



X 2 

V t = XVV rr H V , T > 0, t > (4) 

n 
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Taking the problem — (jSJ) into consideration we have the boundary and the 
initial conditions for equation (0J) in the form 



i dv 

1 n 

dr 



-q Q nt 



t > 



v(r,t = 0)=0, r>0 
To solve the problem (jlj — © we can use the self-similar variables [1-7, 9] 

Br 



v(r,t)=At m f(6), 



e 



(5) 
(6) 

(7) 



where f(6) is a function of 9; A, B, m and p are constants. 
Assuming 



V 



m + 1 



kAB' 



and substituting into equation (jlj) we obtain 

-Ofe -mf = 



//ee + -Je H 



1 



Substituting (jZJ into the boundary condition (jSJ) and taking into account 

m(n +1) (m + 1) 



A- 



n 



n_1 j 2+1 

go-o Tl = A n 

we have the boundary condition for equation Q in the form 

df 



d0 



-1 



(9) 

(10) 
(11) 

(12) 



From equation (JHJ) we find the second boundary condition 

/(0-oo) = O 

Using conditions (jSJ), (jlUjl and (jllj) we get constants A, B, m in form 



.4 



2n „ 



B = h n + 2 q Q n+2 n ™+ 2 
n(2Jfe + 1) 



n + 2 



(13) 



(14) 



Taking expressions (JHJ), (fTUj) and (fTTj) into account we have the boundary value 
problem (0) — Q in the form 



. m + 1 
+ ~fe + — s~ 



1 



6f e -mf = 



f 



" d9 



(15) 
(16) 



/(0-oo) = O (17) 

To find the solution of the problem ((TJ) — (JHJ) we have to solve the problem (fT5J) 
— (117)) . This is the aim of this work. 



2 Method applied. It is known that the velocity of the boundary for the 
nonlinear heat conductivity is finite. Let us assume that 9 = a is the boundary. 
Therefore, at the point a temperature is equal to zero ( f(a) = ), but derivative 
df 



is non-zero 



<:W 



7^0 



We look for approximate solution of the problem (|15J) — ()17)) in the form 



N 



/(*) = £# 



From equation (|TH|) we get 

d 2 f 



d6 



-Pi, 



d9 2 



8=0 



2& 



z!(-l) l A 



Taking /(a) = into account we get @ = 

Substituting (fi"9"j) into equation (JHJ) we have the coefficients 



(19) 



/3i = -an (m + 1) 



1 (m 



n + 1 

1 n(m — 1) (nm + n + 2m) 
12(n + l) 2 a (l + 2n) (m+1) 
1 n(m — 1) (nm + n + 2m) p^ m,?1 ' ) 
48 a 2 ( n + X )3 (i + 2 n ) (m + l) 2 (3 n + 1) 



j(m,n) 



mm 



1) (nm + n + 2m) P t 



(m,n) 



240 a 3 (3 n + 1) (m + If (1 + 2 n) 2 (n + l) 4 (1 + 4 n) 

(m,n) 



where P 4 ' and P 5 ' are polynomials 



(20) 



}(m,n) 



5 nm — n + 7 m — 3 



p (m,n) = 3Q3 + g2 m 2 + 1Q2 ^2 + g 17 ^2 _ 2Q4 ^ 

-238 nm - 70m - 48n 3 m + 12 + 7n 2 - 6n 3 + 31 n 
Using expression (j2"Uj) we can obtain some exact solutions of the problem ()15j) - 

If n + 2m + nm = then (3i = (i > 3). Taking m = consideration 
we have exact solution 



3 1 (m — 1 ) n 

f{9) = -an (m + 1) (a - 6) + - { (a - 

2 4 n + 1 



(21) 



However this exact solution does not satisfy the boundary condition (|16j). This 
solution can be used to solve the Cauchy problem for equation (JTJ). 

If m = 1 then /3j = (z > 2) and we have another exact solution for equation 



f(6) =an (a -6) 
Exact solution satisfies the boundary condition (fTo)). 



(22) 



3 Solutions of the problem (|15|) — (|17|) . Consider the case n > and 
= 1/n, (m = 1). 

f(6) =an(a-6) (23) 
Substituting (J23*j) into the boundary condition (fTTIj) we obtain 

n + l 

= 1 



{n + l)(a 2 nY 



not 



From condition (|2*^ we get the parameter a 



a = [n n (n + 1) 
The values of a and n at k — 1/n are given in table 



(24) 



(25) 



Table 1: 



n 


1 


4/3 


2 


5/2 


3 


4 


a 


1.2599 


1.5299 


2.0598 


2.4586 


2.8619 


3.6840 


n 


9/2 


5 


11/2 


6 


13/2 


7 


a 


4.1025 


4.5256 


4.9528 


5.3838 


5.8184 


6.2561 



Exact solutions of the boundary value problem (|15|) - (fTTj) at A; = 1/n (n > 0) 
are described by formula ([2*3^1. 

Consider the case n = 1 and fc > 0. Approximate solution of the boundary value 
problem (|15|) - ([17]) can be written in the form 



/(f) 

1 (3m+l)(m-l) 



\a (m + 1) (a — 9) + ^(m - 1) (a — 
2 8 



144 



(m + 1) a 



(a-fl) d + 



:j 1 (3m+l)(m-l)(3m-l) / m . 



1152 



(m + 1) a 1 



[a — 



°) (26) 



1 (3m + 1) (m - 1) (201m 2 - 140 m + 11) 5 



172800 (m + l) 3 a 3 

From the boundary condition ()16|) we obtain 



or 



(105147 m 4 + 384822 m 3 + 519188 m 2 + 307082 m + 66161) 



2985984000 ( m + l) 6 (27) 

(9491 + 24237 m 4 + 84342 m 3 + 54602 m + 103808 m 2 ) = 1 
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Figure 1: Comparison of approximate solution (|2*%|) and numerical solution of the 
problem (|TH|) - (fTTj) : 1 — approximate solution; 2 — numerical solution. 

From equation (|27|) we find 



a = 1440 (m + l) 2 / (2548447839 m 8 + 18195239088 m 7 
+55955316456 m 6 + 96920939400 m 5 + 103409323126 m 4 

+69458768096 m 3 + 28562945760 m 2 + 6527038184 m + 627934051) 



Table 2: 



k 





1 


2 


3 


4 


5 


m 


1/3 


1 


5/3 


7/3 


3 


11/3 


a 


1.1762 


0.7937 


0.6222 


0.5211 


0.4532 


0.4039 



The values of a, m and k at n = 1 are given in table EJ Approximate solutions 
of the problem (JT5j) - (fT7|) at n = 1 are expressed by formula (pjj) . 

For the case = approximate solution of the boundary value problem ()15|) - 
(fT7|) takes the form 



2a(a — 



a 



f(0) 



12 



+ 



a — 



a 



144a 



(2f 



0. 



23040a 3 ' 
< ^ < a; 

a < 

where a is a = 1.1762. 

To check approximate solution (|2*Hj) we have compared it with numerical solu- 
tion of the boundary value problem (fTB^) - (fTTj) . The comparison of approximate 
solution (}2*8*j) and numerical one at k = and n — 1 is given on Fig^ Solid line 



5 



is approximate solution and circles correspond to numerical solution of the problem 
(|15|) - (|17p. From FigC]we can see that these solutions are similar. 

Consider the case n = 4/3 and k > 0. Approximate solution of the boundary 
value problem (|T5j) - (fT7j) is given by formula 



f(6) = 2 -a (m + 1) (a - 6) + ~ (m - 1) (a - #) 2 

2 (5m 2 -2-3m)(a-^) 3 
~539 a (m+ 1) 

1 (205 m 3 - 188 m 2 -43 m + 26) (a- 9) A ( 29 ) 
+ 37730 (m+l) 2 a 2 
9 (29055 m 4 + 7078 m + 1199 m 2 - 464 - 36868 m 3 ) (a - 9f 
137997475a 3 (m + l) 3 

Some values of a, m and k at n = 4/3 are given in table El 



Table 3: 



k 





1 


2 


3 


4 


5 


m 


2/5 


6/5 


2 


14/5 


18/5 


22/5 


a 


0.9256 


0.6000 


0.4608 


0.3807 


0.3278 


0.2897 



Consider the case k = and n > 0. From the expression (jT4|) we have 

n 



m 



n + 2 



Approximate solutions of the boundary value problem (|T5j) - (fT7|) at k = 0, n > 
take the form 



L f OL 11? 



1 



n (a — 



+ 



' 2 (n + 1) (n + 2) 
(2n 2 + n-3)n 2 (a - 0) 4 
6«(n + l) 3 (l + 2n) 24 ( n + l) 5 (1 + 2 n) a 2 (3 n + 1) 



n 2 (a-#) J 



1 



1 n 2 (12 + 8n-75n 2 + 12n 5 -77n 3 )(a 



120 



+ . . . 



(30) 



a 3 (n + l) 7 (1 + 2nf (3n + 1) (1 + 4n) 

Using boundary condition (fTB*)) we obtain the value of a. Some values of a and 
n at fc = are given in table 0] 



Table 4: 



n 


1 


4/3 


2 


5/2 


3 


4 


a 


1.4819 


1.1578 


0.7889 


0.6283 


0.5178 


0.3775 


n 


9/2 


5 


11/2 


6 


13/2 


7 


a 


0.3307 


0.2934 


0.2632 


0.2382 


0.2172 


0.1994 
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4 Conclusion. The boundary value problem of the nonlinear heat equation 
for the given flow was considered. This problem was solved using the both numerical 
and analytical approaches. Some exact solutions were found. Approximate solutions 
of the boundary value problem were obtained. Comparison of the numerical and 
the approximate solutions was given. 

This work was supported by the International Science and Technology Center 
(project B1213). 
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